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We present measurements of both /vintegral and /^-differential azimuth quadrupole components 
of two-particle correlations on azimuth (0) and pseudorapidity (77) for unidentified hadrons in 
Au-Au collisions at y/sJm = 62 and 200 GeV. The azimuth quadrupole component is distin- 
guished from 77 -localized same-side correlations by taking advantage of the rj dependence. The 
quadrupole component is related to conventional v>2 measures. Both /^-integral and p t -differential 
results are presented as functions of Au-Au centrality. We observe simple universal energy and 
centrality trends for the /? r integral quadrupole component. p t -differential results are constructed 
using v\ marginal distributions on p,. These results can be transformed to reveal quadrupole 
p t spectra that are nearly independent of centrality. A parametrization of the p t -differential 
quadrupole shows a simple p, dependence that can be factorized from the centrality and colli- 
sion energy dependence above 0.75 GeV/c. 
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1. Introduction 

One of the major physics results of heavy ion collisions at RHIC energies has been the ob- 
servation of a large azimuth anisotropy [p. This phenomenon is typically described in terms 
of hydrodynamics in which there is a pressure-driven expansion of the collision region in non- 
central collisions Hydrodynamics describes both radial flow (a collective transverse velocity 
of particles produced in the collision) and elliptic flow (the azimuth anisotropy produced from the 
asymmetric shape of the original collision system). 

Elliptic flow is typically thought to correspond to the measured quantity V2, defined by the 
second Fourier component of the distribution of particles on azimuth with respect to the angle of 
the reaction place of the colliding nuclei [Q]. This definition works for theoretical calculations in 
which the reaction plane is known. However, in experimental measurements the reaction plane can 
only be estimated by particles in the event, so correlations must be used. 

Two-particle correlations could reliably measure the phenomenon of elliptic flow if there were 
no other physical contributions to the correlation signal [Q]. In practice there are several sources of 
correlations such as (mini)jets, resonances, and HBT effects. There are two major approaches for 
reducing so-called "nonflow" effects. One is to use the pseudorapidity separation between particles 
and the other is to make a many-particle measurement. 

In this analysis we distinguish the azimuth quadrupole from other sources of correlation ge- 
ometrically by using the shape of the correlations on pseudorapidity as well as azimuth. Full 
two-dimensional histograms of the two-particle correlation space are constructed and fitted with 
model functions to distinguish different contributions. 

2. Angular Correlations 

Minimum-bias angular correlations are constructed by considering all possible pairs of parti- 
cles in an event (minus self-pairs). The primary measurement variables of the particles we observe 
are their azimuth angle <p, pseudorapidity 7], and transverse momentum p t . In two-particle correla- 
tions these three variables for each particle define a six-dimensional space (0i,f7ijPfi!02j T 72)/'/2)- 
This space is difficult to work with, so cuts and projections are used to simplify the analysis. 

In this analysis the first cut selects a momentum range from the (pn,Pt2) space. One possible 
choice is to accept the entire momentum range, which yields what we refer to as a /^-integrated cor- 
relation. However, if we want to study the p t dependence of correlations there are several choices, 
since we are dealing with a two-dimensional p t space. The momentum cuts we have chosen are 
marginal distributions in which we restrict the momentum range of one of the particles and allow 
the other particle to have any momentum value. Because of the inherent diagonal symmetry in the 
{ptXiPti) space this produces cross-shaped cuts as seen in Fig. [j] (left panel). 

It is convenient to use transverse rapidity y t instead of transverse momentum p t . y t is a rela- 
tivistic momentum variable defined as 

y, = In {{p t + m,)/m }. (2.1) 

y t is analogous to the usual longitudinal rapidity variable but is defined in the transverse direction 
with y z pa 0. It is useful for studying a large range of momentum values. The definition of y t requires 
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Figure 1: Left Panel: Two example of marginal cuts in y t x y, space. Right Panel: The Fourier decomposi- 
tion of a periodic array of Gaussians as a function of the Gaussian width. 



a mass, but in practice we usually analyze unidentified particles so the pion mass is assumed. When 
used in this fashion y t is a logarithmic transformation of p t that is very close to ln(p f ) but has a 
well-defined zero. 

The correlation measure is constructed from pair densities of sibling pairs (pairs of parti- 
cles taken from the same event) and mixed pairs (pairs of particles taken from similar but dis- 
tinct events). Mixed pairs serve as a reference that can remove detector artifacts in the measure 
Ap / p re f = (p s jb — Pref ) / Pref > w * tn tne mixed-pair reference normalized to the total number of 
the sibling pairs. This is a per-pair measure. At times a per-particle measure such as Ap / y / Pref * s 
more useful since a per-pair measure includes a trivial 1 / n c h trend. 

This approach can be contrasted with the ZYAM (zero yield at minimum) method which is 
applicable only for narrow, well-separated peaks [Q]. The difference is in the normalization of 
mixed pairs to sibling pairs. The ZYAM method imposes the criterion that the minimum of the 
correlation histogram should have value zero and adjusts the normalization in order to achieve that. 
But if the physics signal contains distinct but overlapping peaks then the true correlation value of 
the minimum of the histogram should be nonzero. This normalization also affects the amplitudes 
of the peaks. 

This analysis is based on 14.5 million Au-Au collisions with center-of-mass energy 200 GeV 
and 6.7 million collisions at 62 GeV. Tracks are observed in the STAR TPC with a minimum p t of 
0.15 GeV/c and no upper p t cut. The TPC has 2n coverage in azimuth and this analysis uses tracks 
in — 1 < 17 < 1. The analysis also divides the minimum-bias events into eleven different centrality 
classes as defined in Ref. [B], with nine ~ 10% bins from 100% to 10% and the most-central 10% 
divided into two 5% bins. 

It is useful to compare these correlations to the standard event-plane method for measuring V2 
||], defined by 



, fibs 



v m {EP} 



(2.2) 



(cos(m(y m -Yr)))' 

where V% s = (cos(m(0 — yf m ))), Wr is the true reaction plane angle, \jf m is the estimated event-plane 
angle, and <p is the azimuth angle of a particle in the event. The denominator of this expression is 
known as the event-plane resolution. \jf m is measured as the direction of the corresponding Q vector 



Qm = Y J u( m $i) = Qmu{ m Wm 



(2.3) 
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Here u are unit vectors and Q m is the magnitude of the Q vector. The two-particle correlation 
density V m can be denned by the expression 

V* = t £ u(mfr)u(m$j) = n 2 v 2 m {2}, (2.4) 

where n is the event multiplicity. The second part of the above expression introduces the rela- 
tionship to the two-particle cumulant method. The V m here are strictly Fourier components of the 
total two-particle azimuth density and include "nonflow" effects that might be better described by 
non-Fourier terms. 



If we use Eq. in Eq. 2.4 then we arrive at the following relationship: 



-u^Wm)^- (2-5) 



In Ref. M\ an expression for the event-plane resolution was derived: 



(cos[m(y/- m - Yr)\) ~ y^—^-j^- (2-6) 



Inserting this into Eq. ^5| gives: 

yobs 

V "' ~ n l 7 — = ™m{EP}. (2.7) 

Thus we find that the event-plane method of measuring v%, while often presented using different 
language, is essentially a type of two-particle correlation. 

3. p f -integrated Quadrupole 

p t -integrated quadrupole results have been presented previously in [Q] but are still important 
for understanding the p t - differential results. The systematics they establish provide a basis for fully 
understanding the correlation system. 

The quadrupole component is obtained by doing a free fit of two-dimensional (2D) (t]a ; 0a) 
histograms. It has been found that these histograms can be well-described at all centralities by a 
remarkably simple fit model which works even in proton-proton collisions. The model function 
includes: a same-side 2D Gaussian on (t]a, <Pa), and TjA-independent away-side dipole cos(0a — ft), 
an TjA-independent quadrupole cos(20a)> a ^-independent ID Gaussian on T]a, a narrow same-side 
2D exponential on (t]a 5 0a)> and a constant normalization offset. The model function is expressed 
as 

2 JfOVjLU 2 



-I Ha) ~3 
F = A D cos(<p A -7i)+A Q cos(2<p A )+A e 2 \ a oJ +A x e 



2 ,1 1/2 

+ A 2 e L v " A/ V ' l " Ay J +A 3 . (3.1) 

The quadrupole term is conventionally interpreted as "elliptic flow" [Q]. The other important 
structures in this study are the 2D same-side peak and corresponding away-side ridge, described 
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Figure 2: First Two Panels: p t -integrated quadrupole results in terms of Ap /^/p re f and V2. The solid 
circles and solid triangles are quadrupole results from 2D fitting at 200 and 62 GeV respectively. The open 
circles and open squares are V2{2} and V2{4} results for 200 GeV collisions Jl0[], The upside-down solid 



triangles are 17 GeV results from NA49 [11 1. The open triangles are from projecting the 2D correlations 
to just 0a an d fitting with only Fourier components. The dashed curves are from our parametrization of the 
data. Third Panel: Calculated optical and participant eccentricities as a function of centrality. Right Panel: 
Energy and centrality dependence of the quadrupole using the optical eccentricity for 200 and 62 GeV data 
from this study and 17 GeV results from NA49 [111]]. 



by a 2D Gaussian and the away-side dipole respectively. The same-side 2D peak can contain 
contributions from HBT and resonance decays, but is likely dominated by minijets || ^]. The 
sharp 2D exponential peak mainly describes electron pair production and is not of interest in this 
study. The ID Gaussian on t]a is only present in more-peripheral collisions and is believed to be 
related to participant nucleon fragmentation. Given that it has no 0a dependence it is orthogonal to 
the quadrupole term which is only ^-dependent. 

A more complete description of the away-side structure would be with a periodic array of 
Gaussian peaks centered at it, 3 71, etc. The dipole description applies in the limit where the widths 
of the Gaussian peaks become large. Fig. [I] (right panel) shows the calculated Fourier coefficients 
of such an array of Gaussians as a function of the Gaussian width. In the p t -integrated results the 
dipole appears to be a good description for all centralities. 

The quadrupole component in Ap / y/p Te f is related to the usual V2 measure by [Q] 



Ap[2] 



VPref 



2p v 2 2 (b). 



(3.2) 



Extracted quadrupole parameters are shown in Fig. |2| (first two panels) in terms of both Ap / -y/Pref 



and V2 amplitudes and compared to published STAR V2{2} and V2{4} [10] data and 17 GeV event 



plane results from NA49 flllQ. 

In order to better understand the energy and centrality dependence of the azimuth quadrupole 
it is useful to understand the geometry of the initial collision system which is usually described in 



terms of an eccentricity [ |12| ] 



(yf - (xy 
(y) 2 + {x) 7 



(3.3) 



where x and y are coordinates in the plane perpendicular to the beam axis with x being in the 
reaction plane. Eccentricity is not directly measurable. It is typically estimated based on the 



5 



Azimuth QuadrupoL 



David Kettler for the STAR Collaboration 



centrality of the collision using a Glauber model. The optical model of the eccentricity used here 



is based on a continuous transverse density profile with a Woods-Saxon potential as in Ref. Q13Q. 
In Ref. [0] we show a parametrization of the optical eccentricity for RHIC energies 

, , 1 (3n bin \ - 96 /1136\ 081 

e opt (n bin ) = ^log 10 ^— J -lo^o [— ) • (3-4) 

This is plotted as the solid curve in Fig. |2] (third panel). The dash-dotted curve shows an exam- 
ple of a participant nucleon calculation which is preferred by some to describe conjectured "flow 
fluctuations" [|l4|]. Participant nucleon models give larger eccentricity values at all centralities but 
especially in the very peripheral and central cases. 

In the quadrupole data in Fig. || (first panel) we note that the data at different energies can be 
described by a common shape simply by varying the amplitude. The energy dependence is found 
to be proportional to log(- v /s^v/13 GeV). We then introduce the energy scaling factor 

R(V*m) = log(V%^/13 GeV)/log(200/13). (3.5) 

The complete set of quadrupole data can then be described by 

p {) (b)v 2 {2D}(b) =0.0045R(^J^)e 2 {b)n bin {b). (3.6) 

In Fig. U (fourth panel) we confirm this by observing the linear trend when plotting ( 1 /e 2 ) Ap [2] / y/p re f 
vs. R(y/sm) n bin(b) ■ The parametrization of the quadrupole is very useful because it factorizes the 
energy and centrality dependence. 

We can express the p t -integrated v 2 at centrality b in terms of the p t -differential V2 by 

M b ) = ^j^y J d PtPtPo(Pt,b)v 2 (p t ,b), (3.7) 

where po(p,,b) is the single-particle /? r spectrum [|TJ]. Because the spectrum falls off exponentially 
at larger p t , p t -integrated V2 numbers are heavily weighted toward lower-/?, particles and have 
virtually no sensitivity to V2 above about 0.5 GeV/c. /^-integrated V2 results should not be used to 
conclude anything about the behavior of V2 at higher p t . 



4. /^-differential Quadrupole 

The y t dependence of angular correlations is studied by making cuts on y t and examining 
the 2D histograms in restricted y t intervals. For two-particle correlations it is possible to restrict 
the transverse momentum of each particle independently, so the y t - dependence is inherently two- 
dimensional. While the two-dimensional dependence is interesting in its own right, we construct a 
one-dimensional function on y t to compare to published V2{pt) data. 

We have made nine marginal y t cuts — as described in Sec. || — for each centrality class of 
events for both 62 and 200 GeV collisions. The first bin includes particles from 0.15 GeV — the 
lowest p t we can measure at STAR — to a y t value of 1 .4. Above that there are 7 bins evenly spaced 
in y t with a width of 0.4 units of transverse rapidity. The final bin includes all particles greater than 
y t = 4.2. In p, these bin edges correspond to 0.27 GeV, 0.41 GeV, 0.62 GeV, 0.94 GeV, 1.4 GeV, 
2.1 GeV, 3.1 GeV, and 4.7 GeV. 
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Figure 3: Example of the y t evolution of correlation structures for 62 GeV 40-50% central collisions. The 
plots correspond to y, bins of 1.4 <y t < 1.8, 3.0 < y, < 3.4, and 3.8 < y t < 4.2. 




Figure 4: Example of the y t evolution of correlation structures for 200 GeV 40-50% central collisions. The 
plots correspond to y t bins of 1.4 < y t < 1.8, 3.0 < y t < 3.4, and 3.8 < y t < 4.2. 

In Figs. H and || we show example histograms for 62 GeV and 200 GeV collisions respectively 
at 40-50% centrality. All of the p t -dependent correlations are measured using Ap/p re f which 
simplifies the conversion to V2- 

The fit procedure is similar to that used in the p t -integrated results with one major exception: 
We no longer attempt to model the sharp exponential peak. Instead we simply exclude bins near 
the angular origin from the fit. The excluded bins are at T]a = and A = 0,±jr/12 and at 0a = 
and tja = ±0.08, ±0. 16. This has the effect of making the fits more stable over a wide p t range but 
it can be a problem for the most peripheral bins where the width of the exponential approaches that 
of the same-side 2D Gaussian. 

The major fit parameters are shown in Fig. || for 62 GeV collisions and in Fig. |6| for 200 GeV 
collisions for several centralities. The Gaussian amplitudes and azimuth widths follow expected 
trends for both energies. However, the widths on pseudorapidity seem to be largely independent of 
y t over a large range. 

Quadrupole amplitudes in Ap / p re f can be converted into V2 values by the simple relationship, 
2v\ {2D} = Ap [2] / p re f - However, the marginal distribution produces data of the form 2v2 (p t , b)v2 (b) , 
not 2v\{p t .,b), so the fit amplitude must be divided by /^-integrated V2{b). 

In our fit model the non-quadrupole term with the largest contribution to the second Fourier 
component on azimuth is the same-side 2D Gaussian peak. We can calculate the contribution to 
the second Fourier component of a 2D Gaussian of given amplitude and widths, which we know 
for the same-side peak from our fit parameters. This gives us a direct measure of the so-called 
"nonfiow" contribution to V2, as shown in Fig. [7] for 30-40%, 5-10%, and 0-5% central collisions. 
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Figure 5: Fit parameters for 62 GeV collisions: Same-side peak amplitude, width in eta, width in phi, 
and quadrupole amplitude as a function of y, for 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, and 60-70% 
central collisions. Error bars are for fitting errors only. 




Figure 6: Fit parameters for 200 GeV collisions: Same-side peak amplitude, width in eta, width in phi, 
and quadrupole amplitude as a function of y t for 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, and 60-70% 
central collisions. Error bars are for fitting errors only. 
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Figure 7: Comparison of quadrupole (closed circles) and the second Fourier component of the 2D same- 
side Gaussian (closed upside-down triangles) to V2{EP} results (open circles) [ jl6| ] for 30-40%, 5-10%, 
and 0-5% central collisions. The dark solid curve in the first two panels is the sum of the quadrupole and 
same-side peak terms. 



We see that the 2D Gaussian peak dominates in 0-5% central collisions and that the "nonfiow" (jet 
contribution) is strongly centrality and p t -dependent. 



5. Quadrupole Spectrum 



Consider the case of minimum-bias identified particles results from STAR [17]. In Ref. [18] it 
was shown that by plotting V2 / p t on proper y t for identified hadron species the particles appear to 
come from a single boosted source. This is reproduced in Fig. || (left panel). This is most evident 
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for the more massive particles. Since the pion mass is of the same order as the lowest p t we can 
measure in STAR there are few observed pions in the necessary p t range to see this boost. In Fig. || 
(middle panel) we plot just the Lambdas, which makes the boost more explicit. 

Now consider the centrality-dependent (denoted by b) production of an azimuth quadrupole 
term from a general boosted source. First make the general assumption that the single-particle 
density on y t and (j> can be decomposed into azimuth-dependent and azimuth-independent terms: 

p(y t ,<j>,b)=p (y t ,b)+p 2 (y t ,<j>,b). C 5 - 1 ) 

At this point we do not make any assumptions about the nature of the azimuth dependence. Then 
apply the continuum definition of v 2 : 



^J^ P (y t ^,b)cos[2^-YR)}^Jo n d^p 2 {y,^,b)cos[2^-YR)}_ V 2 

±f^d(t> P (y t ,(t>,b) ^ Po(yt,b) ~Po(yt,b) 



Only the azimuth-dependent term contributes to the numerator, and the denominator is approxi- 



mately the single-particle spectrum [J18p. What we want to study is actually the numerator of this 
expression, which we denote by V 2 . 

To calculate V 2 we need to introduce a boost model. A general boost in nuclear collisions 
should have both monopole (radial flow, Hubble expansion) and quadrupole terms, which is easily 
expressed iny f : 

Ay t (<j>)=Ay t0 + Ay t2 cos(2[<t>-\ir R \), (5.3) 

with Ay t2 < Ay t0 a necessary condition for a positive-definite boost. Using a simple blast-wave 
model with a Maxwell-Boltzmann distribution for a locally-thermalized source the boosted spec- 



trum's azimuth-dependent term then has the form [18]: 



PiiytA) =A 2yt exp{-n 2 [cosh(y t -Ay t (^))) - 1]}, (5.4) 



where p 2 = mo/T 2 . If we insert our boost model into Eq. 1^4] and factor the 0-dependent terms of 



p 2 into the form p 2 (y t ,(j)) = p 2 (y t ) x F^,^) x F 2 (y t ,(j)) where 

Fi(y,,0) = exp{m;[cosh(Ay, 2 cos[20 - y r )}) - l)/T 2 } 

F 2 (y t ,Q) = exp{ iP ;[sinh(Ay r2 cos[2(0 - yr r )])/T z } (5.5) 

then we identify p 2 {y t ) as the azimuth-integrated quadrupole spectrum. Inserting this into the 



numerator of Eq. 5.2 and evaluating the integral yields 



ViMnt^KM. (5.6) 
2l 2 

Finally, we want to isolate p 2 {y t ,b) by taking the unit-integral ratio of measured quantities: 

n , ,x _ V2(y t ,b)/p, _ p 2 (y t ,b) 

Q{yt,b) = „/,>,/,/ r ~ - /,x • (5.7) 
V 2 (b ){\ pt) p 2 (b) 
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Figure 8: Left panel: Identified V2(pt)/pt for pions (closed circles), kaons (open triangles), and protons 
(open circles) for minimum-bias collisions [ |l7| | vs. proper y t for the particle species. Middle panel: A 
magnification of the proton results. Right panel: The quadrupole spectrum for a range of centralities at 
200 GeV. 



The parameters from the blast- wave model present in Eq. 5.6 then drop out in the ratio and Q(y t ,b), 
shown in Fig. || (right panel), directly relates measured parameters to the quadrupole spectrum for 
5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, and 60-70% central 200 GeV collisions. 

We observe approximate centrality-independence of Q(y t ,b) in Fig. || (right panel). We con- 
jecture that there is a centrality-independent Qo(yt) which is well-described by a boosted Levy 
distribution with parameters T 2 = 0.09 GeV, n 2 = 13.8, and Ay t Q = 0.58 (dashed curve), though 
there are significant deviations for the more-central events at low p t . This boost is the same that 
was observed in the case of minimum-bias identified hadrons in Fig. [j§ (left and middle panels). 



6. Quadrupole Parametrization 

A detailed study of the systematics of the p t -dependent quadrupole is possible. First, consider 



rearranging the definition of Q in Eq. 5.7 and use the definition of V 2 from Eq. pJ\to get 



v 2 {2D}(p t ,b) 



p t v 2 {2D}(b) 



Po(b)Qo(p t ) 
Po{Pt,b) 



(6.1) 



The quantity Po{b)Qo(pt) / Po{Pt,b) has a p t dependence described by the ratio of a Levy distribu- 
tion to the single-particle spectrum. This is observed to be exponential for larger values of p t . 
We can then construct a new parametrization of the form 



v 2 {2D\(p,.b) * ( 1 )p ; r : .|2D}-(/>)exp(- P/ /4) < f( Pl .b). 



(6.2) 



where f(p t , b) is a dimensionless factor needed to describe deviations from the exponential form 
at low-/?,. It can be fit to the data with the form 



f(p t ,b) = 1 +C(6)[erf(y, - 1.2) - erf (1.8 - 1.2)], 



(6.3) 



where C(b) = 0.12 — (v — 3.4)/5 — [(v — 3.4) /2] 5 . This parametrization provides a more accurate 
description of the quadrupole term over a wider range of p t and centrality than the Levy distri- 
bution. The factor f(p t ,b) is approximately 1 above about 0.75 GeV/c. Above that point the p t 



10 



Azimuth QuadrupoL 



David Kettler for the STAR Collaboration 



dependence of the quadrupole is entirely described by the factor p t exp(—p t /4). This leads to a 
factorization of the p t and centrality dependence of V2(p t ,b) for these higher p t s. 

This factorization can be combined with the factorization of collision energy and centrality 
dependence of the p t -integrated quadrupole described in Sec. [| for a complete description of the 
azimuth quadrupole component, at least at higher p t . This implies that there is very simple under- 
lying behavior of the azimuth quadrupole. 



7. Conclusions 



Two-particle correlation histograms on azimuth and pseudorapidity have been constructed for 
a wide range of centrality and momentum conditions. These histograms can be fit to study different 
aspects of the physical system. The quadrupole term is closely related to the standard definition of 
V2 but is isolated from r\ -dependent "nonflow" effects. 

Studies of the p t -integrated quadrupole [Q] have revealed simple trends on collision energy 
and centrality. The dependence on collisions energy and centrality can be factorized to produce a 
very accurate description of p t -integrated data. 

We construct p t - dependent histograms using marginal distributions. Published event-plane 



V2(pt,b) data [16] are accurately described by the sum of the quadrupole and 2D Gaussian fit 
components. The /^-dependent quadrupole component is used to construct a boosted quadrupole 
spectrum. The quadrupole spectrum is approximately centrality-independent and is well described 
by a fixed boosted Levy distribution. An accurate parametrization of the p t -differential quadrupole 
exhibits simple scaling above 0.75 GeV/c. 

Quadrupole systematics reveal a system with remarkably simple scaling behavior and a pos- 
sible factorization of the collision energy, centrality, and p t dependence. This seems contrary to 
typical hydrodynamic expectations in nuclear collisions. 
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